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We give a perspective on the Hofstadter butterfly (fractal energy spectrum in magnetic fields), 
which we have shown to arise specifically in three-dimensional(3D) systems in our previous work, (i) 
We first obtain the 'phase diagram' on a parameter space of the transfer energies and the magnetic 
field for the appearance of Hofstadter's butterfly spectrum in anisotropic crystals in 3D. (ii) We show 
that the orientation of the external magnetic held can be arbitrary to have the 3D butterfly, (iii) We 
show that the butterfly is beyond the semiclassical description, (iv) The required magnetic fleld for 
a representative organic metal is estimated to be modest (~ 40 T) if we adopt higher Landau levels 
for the butterfly, (v) We give a simpler way of deriving the topological invariants that represent the 
quantum Hall numbers (i.e., two Hall conductivity in 3D, axy,<Jzx, in units of e^ /h). 



I. INTRODUCTION 

The problem of "commensurability" in almost peri- 
odic systems has been a source of fascination in various 
branches of physics. In this respect, the Bloch states 
in periodic systems in strong magnetic fields have espe- 
cially appealing. Namely, the electron system with a pe- 
riodic potential (or a tight-binding electron system on a 
lattice) in two spatial dimensions(2D) exhibits, when a 
uniform magnetic field is applied, a self-similar energy 
spectrum, known as Hofstadter's butterfly, with fractal 
wave functionad. In this case, the problem is how the 
Bloch electrons react to an external periodicity intro- 
duced by the magnetic field. This problem, raised more 
than twenty years ago by Hofstadter, has been intensively 
investigated. 

The special assets of this problem are: (i) Mathemat- 
ically, Schrddinger's equation reduces to Harper's equa- 
tion with doubly-periodic potential, where one periodic- 
ity is externally governed by the strength of the applied 
magnetic field, B, so that the Bloch states are controlled 
in situ by varying B. 

(ii) In Hofstadter's butterfly (the fractal energy spec- 
trum versus B), each Landau band splits into fractal sub- 
bands, where the integer quantum Hall effect (IQHE) is 
known to occur. The general argument of Laughlina tells 
us that each of the fractally split bands should carry a 
quantized Hall current, while in usual cases each Landau 
band does. Thouless et aa in fact have shown, by actu- 
ally deriving the general formula for the quantized Hall 
conductance in 2D periodic system without any special 
assumption on the periodic potential, that the quantum 



number characterizing the Hall current carried by each 
subband is a topological invariant. 

Now, for both of (i) and (ii) above, two dimension- 
ality is essential, since, first, Schrddinger's equation re- 
duces to doubly-periodic equation only in 2D (while we 
have a triply-periodic potential in 3D) and, second, the 
quantum Hall effect is usually considered to be specific 
to 2D. Nevertheless, we want to ask ourselves here: can 
we have similar phenomena that survive, or even arise 
specifically, in three dimensions (3D)? In the last decade 
several authors have indeed attempted to extend Hof- 
stadter's problem to 3D systemsQa. Usually, however, 
subbands overlap in energy or touch with each other in 
3D, and fractal energy spectra have not been identified. 
On the other hand, HalperinQ has shown, following the 
line of the analysis for 2D, that if the spectrum has an 
energy gap, IQiHE should occur even in 3D. Montambaux 
and KohmotoQ then proposed that 3D Bloch system in 
magnetic fields can realize this situation to derive the 
formula for Hall conductivities. These. results are gener- 
alized by Kohmoto, Halperin and WuQ. 

In our previous letterEfl, we have shown that an ana- 
logue of Hofstadter's butterfly does indeed arise in an 
anisotropic tight-binding lattice {t^ S> ty,tz) when a 
magnetic field is tilted, where the conductive axis is de- 
noted as X. We have calculated the Hall conductivity and 
have shown that each of axy , o'zx , identified to be topolog- 
ical invariants, exhibits an integer quantum Hall effect, 
with a curious one-to-one correspondence with the QHE 
in 2D butterfly. This has been explained with a map- 
ping between 2D and 3D. A crucial point is that the '3D 
butterfly' is specific to three dimensionality rather than 
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a remnant of the 2D butterfly. 

The purpose of the present paper is to address the fol- 
lowing new points: 

(a) to obtain a 'phase diagram' that indicates which re- 
gion in {ty,tz) space should exhibit the butterfly. 

(b) to show that the 3D butterfly is general enough to 
allow arbitrary orientation of the magnetic field. 

(c) to clarify the relation with the semiclassical quanti- 
zation picture in magnetic fields in discussing the details 
of the 2D-3D mapping. 

So the paper is organized as follows. In Scc.||, we re- 
view Harper's equations for the 3D Bloch electrons in 
magnetic fields, and introduce the phase diagram of 3D 
butterfly. In Sec. Ill, we extend the 2D-3D mapping to 
the arbitrary orientation of the magnetic field to discuss 
the influence of ^0 on the butterfly spectrum. In 



Sec. IV, we compare our results with the semiclassical 
approximation. In Sec. VI, we turn to QHE in the 3D 
butterfly, where the Hall conductivity is calculated in a 
much simpler way than in our earlier paper. Experimen- 
tal possibilities to realize the 3D butterfl y is d iscussed in 



DOS 

Sec.0 and the conclusion is given in Sec JVII 



II. THE CONDITION FOR 3D BUTTERFLY 

A. The 2D-3D mapping 

Let us fir-st summarize the results obtained in our pre- 
vious work,EJ where we first suggested that we can have a 
butterfly spectrum in 3D systems. We take a noninteract- 
ing tight-binding electron system in a uniform magnetic 
field, B. Schrodinger's equation is 



(1) 



Here ipi is the wave function at site i, the summation 
is over nearest-neighbor sites, and 9ij — ^ A ■ dl the 
Peierls phase factor arising from the magnetic flux where 
A is the vector potential with W x A = B. 

We consider a 3D simple-cubic lattice in a magnetic 
field B = {0, B sm9, B cos9) assumed to be parallel 
to the yz plane. If we choose the vector potential as 
A = {0,Bxcos6,—Bxsm9), we can substitute in 
eq.(|^) with ipimn = e*"" where l,m,n are the 
site indexes along x, y, z, respectively. Then we have 



tx{Fi-i + Fi+i) - [2tyCOs{2Tr(t>J 
+2tz cos{-2TT(j)yl + iyz)]Fi = 



EF,. 



(2) 



are the transfer integrals between nearest 
neighbors along x,y, z, respectively, and (py or (pz are the 



numbers of fiux quanta penetrating the facet of the unit 
cell normal to y or 2: (Figj 
tentials in (|) V^^^ and 



cell normal to y or 2: (Figj^). Let us call two cosine po- 



If we assume 4>y,4'z *C 1, the periods of these cosine 
potentials are much greater than the original lattice con- 
stant, so that the problem reduces to a continuous ID 
system if we apply the effective-mass approach. One 
can then regard the potential minima of V^(i) as effec- 
tive 'sites', which we call the wells to distinguish from 
the original sites. Each well feels the residual potential 
V^^\ If the well is deep enough, several bound states ap- 
pear at each well, where each state forms a tight-binding 
band, which corresponds to a Landau band. 

Let us define t' as the hopping between adjacent wells, 
where it' gives the width (Harper's broadening) of a Lan- 
dau level. Then eq.(|^) reads, in terms of the wells, 

-t\Jv-i + Jv+i) - 2tzCos[~2Tr{(t>y/<l)z)l' 

+ {<py/(j)z)iyy + Vz\Ji' = EJv, (3) 

where J;/ is the amplitude at the well /'. We can see that 
the equation has the same form as that of the 2D sys- 
tem (eq.(^) with tz =0), which explains why a butterfly 
spectrum can appear also in 3D. If we recall that a 2D 
system exhibits a butterfly when isotropic {tx ~ ty), one 
can predict that the Landau band for which 



t' sa tz 



(4) 



should exhibit a butterfly structure, where the value of 
t' varies from one Landau band to another. 



B. The 'phase diagram' for the 3D butterfly 



We have to satisfy some conditions on [t^ ,ty,tz, (f>y, (j>z) 
to have the 3D butterfly. By systematically examin- 
ing them we can obtain the 'phase diagram', which was 
not obtained in our previous letter. We again assume 
(f>z,4>y ^ Ij and we take tj, = 1 as the unit of energy. 

Let us begin with the pure 2D case {tz = 0), where we 
have only one potential term V^^"^ in Harper's equation. 
If one approximates a well in V^^'^ as a harmonic poten- 
tial, the energy is quantized into Huj = Aircpzy/Uj- The 
condition that there exist more than one levels bound to 
each well (i.e. Landau levels) is then ^Ucu < Aty (depth 
of the well), i.e., 

Tr(t)z < 2^y. (5) 

If the opposite inequality holds, F*^^-' has only a weak 
effect on the energy spectrum and wave functions. 

If we turn tz on, each Landau level feels the residual 
potential V'^'^^ and can split into butterfiy bands. For 
the butterfly gaps to be resolved, however, y^^^ must be 
weak enough so that it does not mix different Landau 
levels. In other words, 4^^, the amplitude of T^'^', must 
be smaller than hio = 47r(/)^^/?r, or 



tz < TTO 



(6) 



We can now construct the phase diagram. Since an 
equivalent discussion can be made when we interchange 
y and z, we must add the conditions 



2 



■K(l)y < 2y/t^ 

ty < n(j)yy^. 



(7) 

(8) 



I = 2tt(I)J, 



The resultant phase diagram is shown in Fig||. The 
phase diagram consists of ID, 2D, 3D butterfly, and 
3D regions. "ID" is the region where neither eq.(|^) 
nor eq.(^ are satisfied. There, both V^^'^ and F^^-* are 
too weak and the energy spectrum can only have Bragg- 
reflection gaps. Since this phase is dominated by t^, we 
call this ID-like. An example, with a Bragg gap, is de- 
picted in Fig.^. "2D-like" is the region where only one 
of eq.(||) or eq.(^ is satisfied. There, one of ty or tz has 
only minor effect and the spectrum is 2D-like in that we 
have only (Harper-broadened) Landau levels. Then we 
come to the "3D butterfly" phases, where eqs.(|[),(^) (or 
eqs.(|^),(||) are satisfied. The remaining region is "3D", 
where all the three hoppings are too large for energy gaps 
to survive. Namely, even when both eq.(^ and eq.(|^) are 
satisfied, the butterfly gaps are wiped out when the other 
conditions are not satisfied, i.e., when the residual poten- 
tial is too large. So the 3D butterfly is allowed to appear 
somewhere in between the 2D and 3D phases. When the 
magnetic fleld is too small, two parabolic phase bound- 
aries coalesce to the axes and the 3D phase becomes dom- 
inant. 

We display typical energy spectrum in each phase in 
Fig.^. In this flgure, panels (b)(c)(d) all belong to the 
3D-butterfly phase, but we can immediately notice that 
the Landau band on which the butterfly appears varies 
according to where we sit on the phase diagram. This is 
because within this phase the 3D butterfly appears when 
t' « tz (eq.^, namely, only the Landau level (not nec- 
essarily the lowest) that has the right width exhibits the 
butterfly. 



C. Estimate of the effective transfer integral 

In our previous letter, we have numerically estimated 
the Landau-level broadening t' for the lowest Landau 
band to predict when 3D butterflies appear. Here we give 
a derivation of the scaling law for t' introduced there. 

Consider first an electron in one dimension in the po- 
tential F(^) = -2ty cos{2Tr(l)J) (no V'-^^) in eq.(|). Here 
we have omitted the phase factor i/y, which has no effect 
on t' . Our aim is to express t' as a function of (t^, ty,(f>z)- 
Harper's equation in this case is 

- i,(Fi_i + Fi+i) ~ 2ty cos{2TT(j)zl)Fi = EFi. (9) 

We have assumed the period of ^^^^ is much greater than 
the lattice constant, so the equation can be cast into a 
differential equation up to a constant, 

- ^ - 2tj, cos(2^0,O^(O - ^^^(0- (10) 



the equation becomes 
d^F . t„ 



dp 



{2t:<Pz)'- 



■coa{l)F{l) = EF{1), 



(11) 



(12) 



where E = E/{2Tr(j)z)^ is the rescaled energy. This equa- 
tion indicates that the form of the wave function depends 
only on ty/(j)j,. 

Since t' depends on the overlap of the neighboring wave 
functions and the amplitude of the periodic potential 
, t' can be written as 



t' = ty f 



(13) 



We again take t^ — 1 as the unit of energy. If we measure 
the length in the well spacing by introducing 



where the function /, which only depends on the Lan- 
dau index, represents the overlap of the neighboring wave 
functions. Hence we have only to calculate t'{ty) for one 
value of (f)z ■ In the limit of large barrier between neigh- 
boring wells, ity/(j>1) — > oo, / vanishes exponentially. In 
the opposite limit of {ty/(j>'j.) 0, wells become shallow, 
so that the transfer to the second- (and more distant) 
neighbor transfers cannot be neglected, and the approx- 
imation should fail. 

This argument can be used for each of the Landau lev- 
els. Higher levels have larger t' since the wave function 
bound in the well is more extended. So, taking account 
of the clear butterfly condition t' « , we find that as tz 
is increased the Landau level in which butterfly appears 
becomes higher, as we have seen in Fig.0. 



III. 3D BUTTERFLY IN AN ARBITRARY 
ORIENTATION OF B 

So far we have assumed that the magnetic field lies 
in the plane (yz) normal to the conductive axis. So the 
next question is whether we can relax this condition for 
realizing the butterfly. About ten years ago, several au- 
thors investigated 3D Hnfetadter problem in arbitrarily 
oriented magnetic fieldslIlEl. They have shown that, even 
if each of three components of magnetic field is not zero, 
Schrodinger's equation can be transformed into a one di- 
mensional tight-binding problem with distant-neighbor 
hoppings. Namely, if we choose a vector potential that 
depends on both y and z, the equation, reduced to a ID 
tight-binding equation, has Uy- and rt^-distant hoppings 
where Uz/uy = 4>z/4'y with Uy, Uz being coprime integers. 
We can adopt this gauge to calculate the 3D butterfiy for 
Bx ^0 here. 

The result for the energy spectra when the magnetic 
field is tilted away from the yz plane {B^ ^ 0) is shown 
in Fig||. Surprisingly, the spectrum is almost identical 
to that for Bx — (Fig.^(a)). In other words, the mag- 
netic field component parallel to the conductive direction, 
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which would promote the motion along that direction, 
has unexpectedly little effect on the 3D butterfly. 

So the puzzle is: why is the 3D butterfly so robust? 
The tight-binding model with distant hoppings cannot 
be reduced to the effective 2D Harper's equation as is 
done here in going from eq.(||) to eq.(||), so we must take 
another approach. 

We consider a 3D butterfly system (t^ ^ ty,tz) in a 
magnetic field {B^, By,B^). We take the vector potential 
A = {BzU — ByZ, BxZ, 0), so that x is a cyclic coordinate 
and the wave function becomes 



H^lrnn — ^ 



(14) 



where {I, m, n) labels {x, y, z). Schrodinger's equation be- 
comes 

''yV^ "-^m — In I ^ ^m+lnj 

mn— 1 Gmn+l) 

- 2tx cos[27r(02m - 4>yn) -f Vx\Gmn = EGmn, (15) 

where 0^ is the number of flux quanta penetrating the 
facet of the unit cell normal to the fi axis. 

We first consider the case Bx = {(j)x = 0). This is 
just the case we considered in the first part of this pa- 
per, but the equation cannot be reduced to a 2D Harper's 
equation due to the different gauge. Nevertheless, we can 
derive an effective 2D Harper equation in the following 
manner. 

In eq.(|l5|), the assumption S> ty,tz means the am- 
plitude of the cosine potential, tx cos[2Tr{cj)zm—(f)yn)+i/x], 
is much greater than the hopping terms along y, z (first 
two terms in eg. (|l5|) ) . So the states are basically bound 
to a 'trough' of the cosine potential. In this sense, the 
state is quasi-one-dimensional. So let us derive an effec- 
tive Schodinger's equation within the trough. 

If we switch off for the moment, the set of equations 
for different n's are decoupled to give 



2tx cos[27r(02m - 0j,n) + i^x]^^^ 



pO/-<0 



(16) 



This is an equation for one variable m \\ y, and the eigen- 
states are bound states in a cosine potential. If tz is 
turned on, the states having different values of n |j z are 
mixed to form a tight-binding band (Fig.^. The band 
width is proportional to tz times the overlap, a, between 
neighboring bound states having a slightly different cen- 
ter positions. The equation can now be reduced to 



atz{Jn^l + Jn+l) 



PO T 



EJn 



(17) 



where J„ is the amplitude of the bound states at n. The 
binding energy in eq.(|l6|) depends on n, since the set 
of values of a steep cosine potential at discrete points 
have a discommensuration that is periodic in n with a 
period 4'z/4'y (i-^-, n ^ n + 4>z/4'y i^i eq.(p^ is equiv- 
alent to TO —s- m — 1). We can check by the numerical 



calculation that the E'^ can be approximated with a co- 
sine function, whose amplitude is denoted by t' . Thus, 
the equation (|l^ becomes 



- atz{Jn-l + Jn+l) 

— 2t' cos[27T{ipy/<j)z)n + const.] J„ = EJn- 



(18) 



This is the effective Harper equation in this approach. 
This equation has the same set of eigenvalues as those of 
eq.(||) as it should. 

If we assume a — l{(t>y <ti (pz), we can prove that t' 



in eqs.(|18D and (|3|) are identical. In fact, eq. (y_8|) has 
the same form as eq.(^ except that t' and tz are ex- 
changed. In the 2D Harper equation, exchanging the two 
hoppings does not change the energy spectrum because 
of the gauge invariance, which concludes the proof. 

The new effective equation can easily be extended to 
the arbitrarily oriented magnetic fields. If Bx is turned 
on, the bound states in a trough feel a magnetic field Bx 
parallel to the yz plane. If the state were purely one- 
dimensional, it would not be affected by magnetic fields. 
So the problem is to what extent we can neglect the ef- 
fect of Bx on the quasi-one-dimensional state. We assume 
0jf <C 02(6* « 0) again. If one considers the state hav- 
ing an energy E above the bottom, the width of the wave 
function can be estimated as ^/E/t^/ipz, while the length 
scale over which the Peierls phase due to Bx changes is 
\/4>x- Thus we can say Bx has a minor effect on the 3D 
butterfly when 



< 



(19) 



In the opposite limit, (py ^ ipziS ~ 90°), we have only 
to replace pz with <f>y in the inequality. In our example 
above, the condition, for the lowest Landau band, be- 
comes (px < ^ziO ~ 0°),(t>x < ipyiO « 90°), which is 



indeed satisfied in Fig.|[ 



IV. RELATION TO THE SEMICLASSICAL 
PICTURE 

The semiclassical picture is often employed in examin- 
ing the motion of electrons in magnetic field. So this sec- 
tion addresses the question of whether or not the present 
approach in terms of Harper's equation can be inter- 
preted in the semiclassical approach. We shall conclude 
that the semiclassical approach cannot describe the but- 
terfly physics. 

In the semiclassical approximation, an electron in k- 
space moves along the intersection of the equienergy sur- 
face with a plane perpendicular to B. If the orbit closes 
in /c-space, we may apply Bohr-Sommerfeld's quantiza- 
tion condition. 



fe_L X dfej 



27r 



(20) 
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which requires that the area enclosed by the orbit must 
be a multiple of 2tt/£'^, where is = \/h/eB is the mag- 
netic length. 

On the other hand, in Harper's equation, the binding 
in wells gives the quantized states. Let us start with iden- 
tifying the correspondence between two quantizations, 
which look different at first sight. In the following, we 



again assume the conditions 



which are our starting point in deriving the effective 
Harper equation. 

The energy dispersion, in the absence of magnetic 
fields, is 

E{k) = —2tx cos k^a — 2ty cos kyb — 2tz cos fc^c, (21) 

where a, 6, c are the lattice constants along x, y, z, respec- 
tively. First we consider a semiclassical orbit around T 
point, around which E(k) can be approximated, up to a 
constant, as 



Eik) = t^{k^af + tyikybf + t,{k,cf 



(22) 



An equienergy surface is then represented by an ellipsoid. 
For B tilted by an angle 9 from z axis, the area enclosed 
by the orbit around F having an energy E is 



S ■ 



ttE 



e^t.:^a'^{tyh'^ +t^c^ tan^ 



cos 



(23) 



Applying the quantization condition, we obtain the quan- 
tized energy. 



En — 



ATTn^t^{ty(t)l + t^(j)l). 



(24) 



Around the X point, (kx, ky, kz) = (0, 0,7r/c) on the 
edge of the Brillouin zone, on the other hand, E{k) can 
be approximated as 



Eik) = tx{kxaf + ty{kyhf - i^vr - fc^c)^ 



(25) 



This time the equienergy surface is hyperbolic along z, 
while elliptic along x^y. The area of an orbit encircling 
X is 



S ■ 



ttE 



(26) 



cos 6\Jtxa?{tyh'^ - tz c?- tan^ &) 
and the quantized energy is 

En - ^T^n^tx{ty<i>l-tz(l>l). (27) 



What corresponds to the quantized energy ( |24| ) or (|27 
in Harper's equation? The potential in Harper's equation 
consists of two periodic terms, 

V^(^^(0 ^ ~2tyCos{2Tr(l)zl + Vy), 
T/(2) (/) = _2i^ cos{-2'K(\)yl + Vz). 

As discussed above, if the condition ty(\)z S> tz'i'y is satis- 
fied, y^^^ is responsible for potential peaks and dips while 



V^^-* provides a slowly varying part. The energy quanta 
are determined by the curvature of the wells. While the 
curvature depends mainly on F^^^ F^^^ also gives an ad- 
ditional contribution. Let us consider the well around a 
minimum of y^^'. We can express the bottom of the well 
in a harmonic form, 

+ ym^t) = ty{2T^^zlf + t,{2T,4>ylf. (28) 

The energy quantum is then given by 



hiOc = ^--K^txa^ityciP;, ^ tz(f>l), 



(29) 



which is identical with eq.(p4[). Around a maximum of 
F^^^ , the potential is written as 

V(^){1) + V'^^\l) = ty(2Tic^zlf - tz{2TTcj)ylf . (30) 

Note that V^^'> and V^"^^ have the opposite curvatures in 
this case. The energy quantum is 



htOr = 47r 



\Jtxa^{ty(f)l - tz(j)l), 



(31) 



which is identical with eq. (p7[) . 

These indicate that a semiclassical orbit in fc-space 
translates to a bound state in a well. In Harper's equa- 
tion, however, a bound state of one well is not an eigen- 
state but hops to adjacent wells, i.e., bound states in dif- 
ferent wells are mixed, which is exactly why each Landau 
level is broadened by the hopping t' . We have seen that 
this effect is essential in explaining the butterfly struc- 
ture. Hence we can conclude that the butterfly is outside 
the semiclassical description. 



V. THE HALL CONDUCTIVITY 

Once the spectrum with a recursive set of energy gaps 
are obtained, a most important question is whether an in- 
teger quantum Hall effect exists on that spectrum. Mon- 
tambaux and Kohmoto have first obtained the quantized 
Hall conductance for three dimensional tight-binding 
electrons in magnetic fields following the line of the re- 
sults for the 2D casdj. In this section, we derive, following 
their arguments, the systematics in the Hall conductivity 
for the 3D butterfly spectrum. Although the results are 
the same as in our earlier paper, we have greatly simpli- 
fied the derivation. 

Let us go back to Harper's equation (H) for the simple 
cubic lattice in rational magnetic fluxes (02;, 0y, 0z) — 
{Q,Py/qy,Pz/qz)- When ipzj'f'y <C 1, we can adopt the 
effective-mass approach for the lattice structure along x, 
and the equation reduces to a differential form, 

Ex{-idx)F[x) - [2tyCos{G:,x + kyh) + 

2tz cos{-GyX + kzc)]F{x) = EF{x), (32) 



5 



where 



r2D = S2Dg2D + t2DP2D- 



(40) 



Gz = 2TT(l)z/a, Gy = 2Tr(j)y/a, 



(33) 



and Ex{kx) is the dispersion along x axis. Now the 
equation is doubly periodic just as in 2D. If we assume 
4>z/4>y — vll (P: 1 :integers), we end up with a Diophan- 
tine equation for the r-th gap, 



(34) 



where s,i are topologically invariant integers, constant 
in each gap. Thus, the number of states per unit volume 
below the rth gap is 



sGz^tGy 27r/6 27r/c 1 _sGz+tGy 
211 /Ly 2-kILzV ~ 



27r6c 



(35) 



where V — LxLyLz is the total volume of tke-system. 

According to the Widom-Stfeda formulalj'tj, the Hall 
conductivity when the Fermi level is in a gap is written 
as 



^ dp 



(36) 



where p is the density of electrons for a fixed chemical 
potential an d Cj^ k is the unit antisymmetric tensor. Sub- 
stituting eq.(p5|) for p, we finally obtain 



eW s t 



(37) 



This is the explicit formula for the Hall conductiv- 
ity derived from Montambaux-Kohmoto's argument. We 
can then calculate the Hall conductivity for an arbitrary 
gap in the 3D butterfly. The result for the Hall conduc- 
tivity is shown in Fig.g(a) for the spectrum displayed in 
Fig.|(a). 

If we compare this with the corresponding plot for 2D 
in Fig.||(b), we can see that there exists a beautiful cor- 
respondence between 2D and 3D cases: azx in 3D corre- 
sponds to (T2D in 2D. The mapping between 2D and 3D 
Harper equations explains this in a transparent manner. 
Let us consider the Fermi level lying in the r'th butterfly 
gap in the A^th Landau level in the 3D butterfly. Since 
one Landau band is composed of q butterfly bands, the 
total gap number r from the bottom is r = Nq + r' . The 
Diophantine equation ( |3^ ) then reads 

r' = (s - N)q + tp. (38) 

On the other hand, according to the argument by Thou- 
less et altl, the Hall conductivity in 2D for the flux 
P2d/'72D is given by 



f2D — — -r^2D, 

h 



(39) 



where ^20 is an integer given by the Diophantine equation 
for the r2Dth gap. 



Recalling the 2D-3D mapping, we can establish a corre- 
spondence 



P 



P2D, 

q < — > 92D, 
r' < — > r2D • 



(41) 



Namely, the correspondence between the topological in- 
variants reads 



N ■ 
t ■ 



S2D, 
t2D- 



(42) 
(43) 



This relation explains the correspondence between azx in 
3D and (T2D in 2D as displayed in the Fig.^. 

We can ask how the Hall conductivities vary when the 
magnetic field is tilted away from the yz plane (with 
Bx ^0). As shown in the Sec. HI, Bx preserves the but- 
terfly structure unless its magnitude is too large. This 
implies that the Hall conductivities, being topological in- 
variants, are not affected by the tilting. 



VI. POSSIBLE EXPERIMENTAL REALIZATIONS 



Let us comment on the possibility of observing the 3D 
butterfly in real materials. In the 2D case, the magnetic 
flux (/) required for the butterfly is the order of unity, that 
is, i? ~ If one assumes an atomic lattice constant 

a — 2 A B ^ 0(10^ )T, which is obviously too huge. In 
the 3D case on the other hand, Fig.|| shows that there 
exists appropriate {ty,tz) for the 3D butterfly no matter 
how small 



')z may be. So, for a given lattice constant, 
the butterfly is more easily realized in 3D. 

In practice, appropriate (ty, tz) become small for small 
values of the fluxes. The energy scale then shrinks when 
the magnetic field is small, so that it will become harder 
to resolve the butterfly structure. For typical quasi- 
ID organic conductors such as (TMTSF)2XEj we have 
tx ■ ty : tz ^ 1 ■ 0-1 : 0.01 with a,b,c ^ lOA, and we can 
estimate the required magnetic field to be 

0.>O.O1, 

which is two orders of magnitude smaller than the 2D 
counterpart, (j)z ^ ^- ipz — 0.01 corresponds to B = 40T, 
which is well within the experimental feasibility. The en- 
ergy scale (the required energji-resolution) is 4:{t' + tz) ^ 
10 meV. In our earlier paperJU we have focused on the 
lowest Landau level in estimating the required field for 
the butterfly, which required a much larger flux (0^ = 
0.1). So the important message here is that, if we employ 
higher Landau levels, the condition is drastically relaxed, 
since the higher levels have larger Harper-broadening. 
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The electron-electron interaction, ignored here, can be 
significant in real materials. Although it may first seem 
that the interaction would smear out the butterfly, we 
can show that it can rather give rise to the butterfly. 
For an organic metal (TMTSF)2X, it has been known 
that the many-body interaction causes a spin density 
wave to become stable in strong magnetic fields, which-is 
called the field-induced spin density wave (FISDW)EJ. 
There, an energy gap opens across the SDW nesting 
wave vector on the Fermi surface, and the hoppings ty 
and give rise to series ot-sja|)S for electron and hole 
bands around the SDW gapllill. In (TMTSF)2X where 
tx : ty : tz ^ 1 : 0.1 : 0.01, happens to be too small 
to play a role, so that the spectrum is 2D-like, i.e., has 
only Landau levels around the main gap. We can then 
propose that, if there are quasi-lD materials that have 
ty « tz (rather than ty 3> tz), the spectrum is expected 
to have 3D butterfly gaps according to the present theory. 
This mechanism is interesting, since the butterfly forma- 
tion, being SDW related, is around Ep by construction. 
The detail will be published elsewhere. 



VII. CONCLUSIONS 

We have investigated the energy spectra and quantum 
Hall effects in anisotropic three-dimensional Bloch elec- 
trons in magnetic fields. We have (i) obtained a phase 
diagram versus (ty ,tz) to identify the region in which 
the butterfly appears, (ii) We have extended our the- 
ory to the arbitrary orientation of B, and proved that 
a change in the orientation has surprisingly little effect 
on the 3D butterfly, (iii) We have discussed the relation 
between our method and the semiclassical approach, and 
found that while there is a clear correspondence in Lan- 
dau's quantization, the butterfly structure is beyond the 
semiclassical picture, (iv) The required magnetic field 
for a representative organic metal (TMTSF compound) 
is estimated to be modest 40 T) if we adopt higher 
Landau levels for the butterfiy. (v) We have finally given 
a simpler way of deriving the topological invariants that 
represent the quantum Hall numbers (i.e., two Hall con- 
ductivity in 3D, cTxj/, cfzx, in units of jh). 

M.K. would like to thank the JSPS Research Fellow- 
ships for Young Scientists for financial support. 

FIG. 1. The geometry in a 3D lattice in a tilted magnetic 
field. 

FIG. 2. The phase diagram that indicates the char- 
acter of the energy spectrum in anisotropic 3D systems 
{ty,tz <S ta; = 1) against ty and for the given values of 
the magnetic fluxes (pyj^iz as deflned in Fig.nl. The 3D but- 
terfly region is shaded. 



FIG. 3. Typical energy spectra plotted against tan6' for 
various positions in the phase diagram (indicated on Figji] as 
(a)-(e)). Here we have fixed {ty = 0.5, (j)^ = 1/10) to vary 
t./4>l. We can see that the butterfiy region shifts upward as 
tz is increased, which confirms the above observation since the 
width of the Landau band is larger for higher Landau bands. 

FIG. 4. Typical ID-like energy spectrum for 

tx '■ ty : tz — 1 : 0.01 : 0.01 plotted against 9 in magnetic 
fields {(jjy, (f)z) = l/5(sin 9, cos 9). The gaps are due to Bragg's 
refiection by V^^' and V^^'. 

FIG. 5. The energy spectra of the 3D system 

tx ■ ty : tz = 1 : 0.1 : 0.1 plotted against the tilting an- 
gle in a magnetic field y^l + ^l = 1/5 for various tilting 
orientations of B (as depicted in the inset). 



FIG. 6. The potential and wave functions along y are 
schematically depicted for eq.(|l7|). The grayscale represents 
the value of the potential. The eigenstates bound at different 
positions in the trough are mixed by tz- 

FIG. 7. Semiclassical electron orbits in fc-space are 
schematically shown. We have assumed tx > ty > tz. 

FIG. 8. (a) Hall conductivities {(Jxy,cTzx) ~ —{e^ /h){s,t) 
are plotted on a 3D butterfiy, where we display the topological 
invariants {s,t) for each gap. (b) The corresponding plot for 
CT2D = — (e^//i)t2D on the 2D butterfiy. The region enclosed 
by a dashed line in (a) corresponds to the 2D butterfiy, where 
t in (a) corresponds to t2D in (b), while s in (a) to S2d in 
Eq.(g|) 
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